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The Lyapunov spectrum corresponding to a periodic orbit for a two 
dimensional many particle system with hard core interactions is dis- 
cussed. Noting that the matrix to describe the tangent space dynam- 
ics has the block cyclic structure, the calculation of the Lyapunov 
spectrum is attributed to the eigenvalue problem of 16 x 16 reduced 
matrices regardless of the number of particles. We show that there 
is the thermodynamic limit of the Lyapunov spectrum in this peri- 
odic orbit. The Lyapunov spectrum has a step structure, which is 
explained by using symmetries of the reduced matrices. 



1 Introduction 

Just as low dimensional chaos is revealed by a single positive Lyapunov exponent 
which indicates exponential sensitivity to initial conditions, systems with phase 
space of very high dimension can be characterized by their Lyapunov spectra 
which give information about many possible instabilities in the system. As 
a concrete example, we consider the system consisting of N disks with hard 
core interactions and periodic boundary conditions. This is a surprisingly good 
model of a fluid [Q, and yet is sufficiently simple that ergodic properties may 
be established under fairly general conditions 0] . 

The following is a brief review of the study of Lyapunov spectra in such 
systems; for more details and references in numerical work, see Ref. || and in 
analytical work, see Ref. 0. It has been known for some time that the Lya- 
punov exponents of Hamiltonian systems come in plus/minus pairs, that is, the 
spectrum is symmetric about zero Non-equilibrium extensions were shown 
to exhibit symmetry about a point other than zero R), leading to the dis- 
covery that these extensions also contain hidden Hamiltonian structure ||, ^j. 
Also known for more than twenty years is the algorithm for numerical compu- 
tation of Lyapunov exponents due to Benettin and others |l0| , O], Later, a 
constraint method was introduced More recently, the effects of the hard 

collisions have been properly taken into account [|[ [l| . The existence of a ther- 
modynamic limit in Lyapunov spectra, that is, that the spectrum retains its 
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shape as the number of particles increases, has been put forward using random 
matrix approximations |I3, numerical evidence |l5| and mathematical argu- 
ments Jl6| , but recent numerical work has suggested in contrast, a logarithmic 
singularity of the largest Lyapunov exponent with the number of particles |l7j . 
Lyapunov spectra for diatomic molecules (represented by dumb-bells) show an 
explicit separation of the rotational and translational degrees of freedom if the 
departure from sphericity is small enough p8| . Finally, and of particular in- 
terest to us, careful simulations of sufficiently large systems have revealed a 
step structure in the Lyapunov spectrum for the smallest positive Lyapunov 
exponents H, OL |l8). The above references give an incomplete description in 
terms of "Posch Lyapunov modes" , phase space perturbations corresponding to 
these small Lyapunov exponents which are approximately sinusoidal in position 
space. Recently Eckmann and Gat have suggested an explanation of the Lya- 
punov modes of a one dimensional system using a random matrix approximation 
of the Lyapunov spectrum [jl9| . In this paper we study step structure in the 
Lyapunov spectrum of a two dimensional many particle system without making 
any approximations, however we are restricted to periodic orbits. 

Periodic orbit theory ^lj has proven very useful for investigations of the 
corresponding low dimensional system, the periodic Lorentz gas |22|, |24|, , 
computing properties such as the diffusion coefficient. These methods cannot 
generally be applied directly to high dimensional systems due to the difficulty 
of finding all the periodic orbits, although a notable exception is the Kuramoto- 
Shivashinsky PDE in a regime where the effective number of degrees of freedom 
is small Q. However, periodic orbit arguments have been used to justify 
thermodynamic results such as non-negativity of the entropy production p7{ 
and the Onsager relations [^8| without explicitly finding any periodic orbits. 
Periodic orbits of spatially extended systems in the form of coupled map lattices 
have been considered previously |2^, [30| , leading to block cyclic matrices similar 
to those observed in this paper. 

Here we apply methods similar to those to Ref.|50| to reduce the problem 
of computing the Lyapunov spectrum to that of finding the eigenvalues of a 
relatively small (16 x 16) matrix. A step structure is observed, which is related 
to the symmetries of the system. Our formalism shows that, at least in the case 
of this periodic orbit, the thermodynamic limit of the Lyapunov spectrum holds 
exactly 



2 Lyapunov Exponents of Periodic Orbits for 
the Many-Particle System 

The system which we consider in this paper is a two-dimensional Hamiltonian 
system consisting of N disks, interacting only by hard core collisions. We choose 
units such that the mass and radius of the particles are one. We write the 
position and the momentum of the j-th particle as and p ■ , respectively, and 
for a later convenience we represent the phase space vector r as a column vector 
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((ZiiPii " " > <1niPn) T witli T the transpose operation. 

The dynamics of such a many-particle system is simply separated into the 
free flight part and the collision part, and the tangent vector Sr n of the phase 
space just after the n-th collision occurs is related to the tangent vector STq at 
the initial time by 5T n = M n 5rQ with the AN x AN matrix M n represented as 



M n = M^M^M^M^ ■ ■ ■ M[ c) M[ f) . 



(1) 



( f) 

Here M- is the AN x AN matrix to specify the j-th free flight dynamics, and 
is given by 



M 



(/) 



-(!)> 
J N I 



(2) 



where Diag(Xi, X2, • ■ • ,X{) means the matrix on whose diagonal are the matrix 
blocks Xi,X%, ■ ■ ■ 1 Xi for an integer I, and Up is defined by 



-d) 



I 



(3) 



where 7 and are 2x2 identical and null matrices, respectively, and Tj is its 
corresponding free flight time. On the other hand, is the AN x AN matrix 
to specify the j-th collision of particles. For the simplest case in which the j-th 

(c) 

collision involves only the kj-th and the lj-ih particles colliding, the matrix 
is given as the block matrix 



/ / 



M 



G, 



i J 



(4) 



where I is the 4x4 identity matrix and is put in the part of ' • of the above 
representation, and Fj (Gj) is the (kj, kj) and (lj, lj) ((kj, lj) and (lj, kj)) block 
matrix elements, and the 4x4 null matrices are put in the other elements. Here 
4x4 matrices Fj and Gj are defined by 




(5) 
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Gj = J 3) 7 2) (6) 



(2) (3) 

with the 2x2 matrices L - } and L K - ' defined by 

Lf s ni nJ (7) 



n 



(Note that all vectors in this paper are introduced as column vectors, so for 
example, njApW is a scalar and rijAp^ T is a matrix.) where rij is the unit 
vector pointing from the center of the kj-th disk to the center of the lj-th disk 
at the j-th collision, and Ap^' is the momentum difference p l} —p kj just before 
the j-th collision. 

Now we consider the case that the movement of the system is periodic in 
time, so that the condition r n = Po is satisfied for a integer n p . In this case 
the Lyapunov exponents Xj, j = 1, 2, • • • , AN of the periodic orbit are defined 
by Xj — (l/t(n p )) In \nij(n p )\ in terms of the absolute value of the (generally 
complex) eigenvalue mj(n p ) of the matrix M„ p , where t(n p ) = Y^j=i T j ^ s ^he 
period of this orbit. 

We put the set of the Lyapunov exponents in descending order, namely 
Xi > A2 > ■ • • > A47V • It is well known that in the time-independent Hamiltonian 
system the Lyapunov exponents satisfy the pairing rule |^| , namely the condition 
Ai + X^n = A2 + A4JV-1 = • • • = X2N + X2N+1 — 0. Noting this fact, thereafter 
we consider only the first half Ai, A2, • • • , X2N of the full Lyapunov exponents, 
and refer their set as the Lyapunov spectrum in this paper. 



3 Periodic Orbit Model and the Reduced Ma- 
trix 

By using the method given in the previous section, we calculate the Lyapunov 
spectrum for the periodic orbit illustrated in Fig. ^. In this periodic orbit each 
particle moves in a square orbit with the same absolute value of momentum 
and the same direction of rotation and with a constant free flight time r (So 
the period of the orbit is t(4) = At.). We put 2N\ (2N2) as the number of the 
particles in each horizontal line (each vertical line) so that N = 4NiN2- We 
impose periodic boundary conditions, thus requiring the number of particles in 
each direction to be even. 

In this periodic orbit model, there are only two types of collisions illustrated 
in Fig. |^, in which the trajectories of two colliding particles are drawn with 
their moving directions shown by the arrows. Now we consider the dynam- 
ics of particles for a free flight plus one of these two collisions. Such a dy- 
namics is described by the matrix multiplications F\ = FiDiag(L^ 1 \ L^) and 
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Figure 1: Periodic orbit of the many particle system. The circular dots show the 
positions of particles at the initial time, and the solid lines give the subsequent 
path in the direction of the arrows. 



GisGiDiagCL^.L^) (F 2 



F 2 Diag(4 1} 



4 X) ) and G 2 = G 2 Diag(L 2 ±J , L^>)) 
' By using Eqs. (§),(§), ® and 



(i) r (lh 



corresponding to the left (right) orbit in Fig. 
(Bp these matrices are simply given by 
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Figure 2: Two types of particle collisions. The solid lines give the path of the 
particles, which move in the direction shown by the arrows. 



where one of the components of the collision vector is zero, namely n = (1,0) T 
or (0, 1) T , and v/V2 is the speed of the particles. 

The AN x AN matrices j = 1, 2, 3, 4 are represented as 



M[ c) M[ f) = Diag(Pi, Pi, P u P 2 , ■ ■ ■ , Pi, P 2 ) 



(13) 



M 2 (c) M 2 (/) 



/ Qo Qi 

Qi Qa Q2 

Q2 Qa Qi 

Qi Qo Q2 



Q2 \ 



V Q2 



Q2 Qo Qi 
Qi Qo ) 



(14) 



M^M { J ] = Diag(P 2 , P l5 P 2 , Pr, ■ • • , P 2 , P) 



(15) 



/ Qo Q2 
Q2 Qo Qi 



M[ c) M[ f) 



Qi \ 



!l WO W2 

Q2 Qo Qi 



V Qi 



Qi Qo Q2 
Q2 Qo J 



(16) 



where the 8N1 x 8A^i matrices Pj and Qj are defined by 
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( h Gi 

G\ F 



Pi 



Fi 



Gi 

h 



Fi 



Fx 



Fx Gx 
Gx Fx J 



( Fx 



Pi = 



Fa Gx 
Gx h 



Fx Gx 
Gx Fx 



Gx \ 



V Gx 



Fx Gx 
Gx h 



Fi J 



Qo = Diag(F 2 , F 2 , F 2 , F 2 , ■ ■ ■ , F 2 ,F 2 ) (19) 
Qx =Diag(G 2 ,0,G 2 ,0,---,G 2 ,0). (20) 

Q 2 = Diag(0,G 2 ,0,G 2 ,---,0,G 2 ) (21) 

In order to get these expression of the matrices the first horizontal row of par- 
ticles is numbered 1,2, ■ ■ ■ ,2Nx from left to right, the second row numbered 
2Nx + 1, 2Nx + 2, • • • , 4Nx and so on until the last row numbered (2N 2 - l)2Nx + 
l,(2N 2 -l)2Nx+2,---,4N 2 Nx. 

The matrix M Hp , whose eigenvalues lead to the Lyapunov spectrum, is given 

by M np = M^M^M^M^M^M^M^M^. 

It is very important to note that this periodic model is invariant with re- 
spect to translations horizontally or vertically by the distance corresponding to 
two particles. This feature is reflected as block-cyclic structures in the matrix 
M n . By using this fact, as shown in Appendix A, we can show that the eigen- 
values of the matrix M„ p are equal to the eigenvalues of the 16 x 16 matrices 
M(2irnx/Nx,2irn 2 /N 2 ), n x = l,2,---,Nx, n 2 = 1,2,---,N 2 with the matrix 
M(k, I) defined by 
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M(k,l) = 

/ Si(-k,-l) S 2 (-k,-l) Tx{k,l)e- il T 2 (k,l)e^ k +^ \ 

S 2 (k,l) Si(k,l) T 2 (-k,-l)e ik T x (-*,-!) 

Tx(-k,-l)e il T 2 (-k,-l)e a Sx(k,l) S 2 (k,l) e -* k 

\ T 2 (k,l) T x (fc,Z) S 2 (-k,-l)e ik Sx(-k,-l) J 

Here, Si(fe,Z), «5 2 (fc,Z), Ti(fc,Z) and T 2 (fc,Z) are denned by 



Sx(k,l) = (F 2 A) 2 + (G 2 Gi) 2 + (F 2 G!) 2 e 4fc + (G 2 Fi) 2 e il (23) 
S 2 (M) = F 2 F\F 2 Gx + G 2 G\G 2 F\ + F 2 G\F 2 F\e ik + G 2 FxG 2 Gxe il (24) 
Ti(fc, ee F 2 GiG 2 G x + G 2 FxF 2 Fx + G 2 GxF 2 G\e lk + F 2 F 1 G 2 F 1 e il (25) 

T 2 (k, I) ee F 2 G x G 2 Fx + G 2 F x F 2 Gx + G 2 GxF 2 F x e lk + F 2 F 1 G 2 G 1 e a . (26) 

In the next section we investigate the Lyapunov spectrum of this periodic orbit 
model by using the eigenvalues of the matrix A4(k, I). 

4 Lyapunov Spectrum and its Step Structure 

The reduced matrix ( p2| ) is useful, not only to reduce the calculation time to get 
the full Lyapunov spectrum, but also to allow us to consider some properties of 
the Lyapunov spectrum itself. One of important results obtained by such a con- 
sideration using the reduced matrix is the existence of the thermodynamic limit 
in the Lyapunov spectrum. It should be noted that the NxN 2 matrices Ai(k, I), 
with k = 2-kux/Nx, I = 2ttti 2 /N 2 n x = l,2,---,Nx and n 2 = 1, 2, • • • , N 2 have 
the same form of the matrix given by Eq. (p2|), so in the limit Nx — * +oo 
and N 2 — > +oo the Lyapunov spectrum is given through the eigenvalues of the 
matrices A4(k, I), k £ [0, 27r) and / £ [0, 2tt). This also implies that the max- 
imum Lyapunov exponent takes a finite value in the thermodynamic limit in 
this periodic orbit model. 

Now we calculate the Lyapunov spectrum in our model. Fig. || is the Lya- 
punov spectrum in the case of v = 1.8, r = 2.3, Nx = 11 and N 2 — 9 (So the 
total number of particles is N = AN\N 2 = 396). One of the remarkable features 
of this Lyapunov spectrum is its step structure. It is important to note that 
some steps of the Lyapunov spectrum are explained by using symmetries of the 
reduced matrix (|2^). Actually we can show (at least numerically) that in the 
case of Fig. ||the Lyapunov exponents calculated using the matrix A4(k, I) are 
invariant under the transformations k — * 2ir — k and Z — > 27r — Z, leading to steps 
in the Lyapunov spectrum. Performing both of these transformations at once 
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Figure 3 

Figure 3: Lyapunov spectrum in the case of {Ni,N 2 ) — (11,9). (a) Full scale, 
(b) Large j part. 



has the effect of taking the complex conjugate of M. so the result is obvious, 
however it is not immediately obvious that the spectrum should be invariant if 
k and I are transformed separately. 

We can investigate the relation between a symmetry of the system and the 
step structure in the Lyapunov spectrum another way. For this purpose we 
consider the Lyapunov spectra in the square case and the rectangular case with 
the same number of particles. Here, the square system has the symmetry for 
the exchange of the vertical and the horizontal directions, which the rectan- 
gular system does not have. Fig. ^ is the Lyapunov spectra in the case of 
(Ni,N 2 ) = (12,12) (the upper two graphs) and in the case (Ni,N 2 ) = (24,6) 
(the lower two graphs) with v = 1.8 and r = 2.3. These graphs show that there 
is not a remarkable difference in the global shapes of the Lyapunov spectrum 
in these two cases, but the square system has (even twice) longer steps in the 
Lyapunov spectrum than in the rectangular system, shown in Figs. || (b) and 
(d). These longer steps come from an additional symmetry in the square system. 
Actually we can check numerically that the Lyapunov exponents obtained from 
the reduced matrix ( |22] ) are invariant under k <-> I, which leads to more degen- 
eracy in the case Ni = N 2 = 12 than when (Ni,N 2 ) = (24,6). For Ni = N 2 
all that is required is that n\ and n 2 are interchanged. For (Ni,N 2 ) = (24,6) 
degeneracy by this mechanism only occurs if n\ is divisible by 4 which is much 
rarer. 

It should be noted that some steps of the Lyapunov spectra are too close to 
be distinguished in Fig. ||. For example, if we can investigate more precisely, 
we can see 5 different steps in the long flat part of the Lyapunov exponents A., 
in j e [1018, 1097] (in j £ [1086, 1125]) in the Lyapunov spectrum in the case 
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Figure 4 



Figure 4: Lyapunov spectra in the cases of a square system (Ni,N 2 ) = (12, 12) 
((a) Full scale, (b) Large j part.) and a rectangular system (Ni,N 2 ) — (24,6) 
((c) Full scale, (d) Large j part.). 



of (7Vi,7V 2 ) = (12,12) (in the case of (N l7 N 2 ) = (24,6)). Similarly, most of 
the Lyapunov exponents Xj, j <G [1106,1152] (j £ [1130,1152])) in the case of 
(Ni,N 2 ) = (12,12) (in the case of (N l7 N 2 ) = (24,6)) are not zero, just too 
small for the scale of the plot. 



5 Conclusion and Remarks 

In this paper we have investigated the Lyapunov spectrum of a periodic orbit of 
a two-dimensional system, which consists of many disks with hard-core interac- 
tions. The system has a rectangular shape (or a square shape in a special case), 
and we used periodic boundary conditions. By the block cyclic structure of 
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the matrix, which describes the tangent space dynamics, the calculation of the 
Lyapunov exponents is simplified to the eigenvalue problems of reduced 16 x 16 
matrices however many particles. The reduced matrix is used to consider the 
relation between such a step structure and symmetries of the system, and to 
show the existence of the thermodynamic limit in the Lyapunov spectrum. In 
particular we showed that the difference of the aspect ratio of the system ap- 
pears in the stepwise structure of the Lyapunov spectrum rather than in the 
global shape of the Lyapunov spectrum. 

The Lyapunov spectrum of the periodic orbit discussed in this paper shows 
a step structure, but it should be noted that this step structure is different 
from the step structure obtained by the numerical work [0. More systematic 
investigations of the Lyapunov spectra of periodic orbits may be requested to 
explain the numerical results. 

It should be emphasized that there are many periodic orbits in which we 
can calculate their Lyapunov spectra in many-hard-disk systems by using the 
block cyclic technique shown in this paper. On the other hand we are still far 
from the position where we can calculate the general Lyapunov spectrum of 
the many hard disk system by using the periodic orbit expansion technique. 
One of the problems is that in many particle system we could not know how to 
systematically find all of the periodic orbits whose periods are less than a given 
length. In addition, the periodic orbits of many-panicle system are typically 
distributed continuously, not isolated from each other. These problems make 
the application of the periodic orbit expansion technique to the many-particle 
systems difficult, and not yet solved. 
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A Block Cyclic Structure and the Reduced Ma- 
trix 

In this appendix we show that the eigenvalues of the matrix M n are equal 
to the eigenvalues of the matrices A4(2Trni/Ni,2Trn,2/N2), n\ — 1,2, •••,iV 1 , 
n 2 = 1,2,-..,JV 2 , given by Eq. ©. 

First we calculate the matrix M„ p . We multiply the matrices 
j = 1, 2, 3, 4 given by Eqs. (|f), @, @ and (H), and obtain 
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M nv 



( K K x 
K 2 K K\ 
K 2 K 



where Kj, j = 0, 1, 2 are defined by 



K 2 \ 



K 2 K K x 
K 2 K ) 



(27) 



K 



L02L01 + L 2 \L\\ + L\\L 2 i 
LqiLu + L 22 Lqi 



Lq 2 Li 2 + L 2 iLq 2 



Lm L 



01-^02 



L27L1 



(28) 



K x = 



r 2 







+ ^12-^01 -^12 



(29) 



1 L lx LhLq 2 + Lq 2 L 22 
* 2 =< 6 L\ 2 



(30) 



with the matrices Ljk = QjPk and %N\ x 8-/V1 null matrix 0. Eq. (g7|) implies 
that the matrix M n has the block cyclic structure, so it is block-diagonalized 
by the orthogonal matrix U(16N\, N 2 ) introduced through 



/ j.g27rilxl/fc j.g27ril 



U(j,k) 



x2/fc 



r „2iri2xl/fc j 27ri2x2/fc 



j. e 2TTilxk/k \ 
J.g2iri2xfc/fc 



(31) 



1 t 2-nikxllk j 2irikx2/k j 2-rrikxk/k , 



with the j x j identical matrix Ij so that we obtain the matrix U(16Ni, N 2 y 
M np U(16N u N 2 ) = Diag(A(27rl/iV2), A(2tt2/N 2 ) 1 A(2%N 2 /N 2 ) with j mean- 
ing to take its Hcrmitian conjugate. Here the matrix A(l) is defined by 



A(l) = Ko + K ie u +K 2 e- 



(32) 



The eigenvalues of the matrix M n are equal to the eigenvalues of the matrices 
A{2im 2 /N 2 ), n 2 = l,2,---,N 2 . 

The matrix A{1) is represented as 



A(l) = 



(33) 



12 



where B^\l) is defined by 



BV>(1) 



( B<j\l) B[ j \l) 

BP(l) BP® B[ j \l) 

B«\l) B<j\l) B[ j \l) 



V B[ j \l) 



B«\l) B<j\l) B[ j \l) 

B%\l) Bjf\l) ) 



Here the matrices B^ (I), j = 1, 2, 3, 4 and k = 0, 1, 2 are defined by 



F2F1F2G1 - 



(G 2 Gi) 2 + (G 2 F 1 ) 2 e~ 



G2G1G2F1 + G2F\G2G\e 1 ' 



F2F1F2G1 + G2G1G2F1 + G2F\G2G\e 
(F 2 F 1 ) 2 + (G2G0 2 + (G 2 #i) 2 e a 



B«(i) 











F2GIF2FI (F 2 Gl)' i 



R (i)m_ f (^GO 2 F 2 G 1 F 2 F 1 

b 2 W = I 6 6 



F 2 F l G 2 F l + {G 2 F 1 F 2 F 1 + F2G l G 2 G l )e' a 
G 2 GiF 2 F\ 

GjGiF 2 Fie- a _ 
F2G1G2G1 + G2F1F2F1 + F2F\G2F\e 



B (2) m= / GsGi^Gie- 4 ' 6 

1 U \ F 2 G 1 G 2 F 1 + G 2 F 1 F 2 G 1 + F 2 F 1 G 2 G 1 e- u 6 



s ( 2)(z) ^ / F 2 F 1 G 2 G 1 + (Ga/ifbGiJ- F 2 GiG 2 F 1 )e~ il 

V g 2 G\F 2 G\ 



(3), n f F 2 F l G 2 F 1 + (F 2 G 1 G 2 G 1 + G 2 F 1 F2p 1 )e il 

U - I, F 2 G 1 G 2 F 1 + G 2 F x F 2 Gx + F 2 F 1 G 2 G 1 e il 

F 2 F\G 2 Gi + (F2G1G2F1 + G 2 F\F 2 G{)e 
F 2 G\G 2 Gi + G 2 F\F 2 Fi + F 2 F\G 2 F\e i 
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B^(l)= ( ^zGi-F^Gie 1 ' G 2 G\F 2 F\e % j , ( ., 



y F2G1F2F1 

F2G1F2F1 
(F 2 F 1 ) 2 + (G 2 G 1 ) 2 + (G 2 F 1 ) 2 e- il 



(44) 



B W(I) s f - P^i) 2 6 ] (45) 

V F 2 F\F 2 G X + G2G1G2F1 + G 2 FiG2Gie 1 



r( 4 ) ('/">= ( F 2 F-yF 2 G x + G 2 G\G 2 F\ + G 2 F-yG 2 G x e % \ 
2 1 j " I 6 (F 2 Gi) 2 1 



with the 4 x 4 null matrix 0. It follows from Eqs. fl3|), (|3j) and (|34j) that 
Diag(C/(8, ATi), 17(8, ATi))U(Z) Diag(f7(8, N t ), U(fi, N,)) 



(7(3) (;) crW(0 



(47) 



where CW(Z), 3 = 1,2,3,4 are defined by C^(l) = Diag( Z>W(27rl/JVi, Z), 
DW(2ir2/Nx,l), • • • , DU)(2nNx/Ni,l) ) with £>W(jfe, Z) ee +B^V fc -t-B^e-^. 
The matrix .M(Zc, Z) is introduced as 

Af(fcZ)-^ (1)(fe ' Z) (48) 



which is equal to Eq. (|22j). Therefore the eigenvalues of the matrix A{1) are 
equal to the eigenvalues of the matrices M.(2Trni/Ni,l), n x — 1, 2, • • • , N\. It 
implies that the eigenvalues of the matrix M Up are equal to the eigenvalues of 
the matrices M{2-Kni/N u 2m 2 /N 2 ), n x = 1,2,-- -,N X , n 2 =1,2,---,N 2 . 
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